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Abstract 



of NN scattering must be incorporated in order for one to obtain satisfactory convergence for chiral effective 

CN ■ field theory. Without the proposition of perturbative one-pion exchange, an effective field theory with the 

spin-0, isospin-1 dibaryon is developed to account for this fine-tuning, and is demonstrated up to 0{Q^) 

where the leading irreducible two-pion exchange arises. The approach shown in the paper results in rapid 

O ■ convergence of the ^So partial- wave amplitude. 
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INTRODUCTION 



In the framework of chiral effective field theory (ChEFT), Weinberg's original power counting 
(WPC) [jtSI] for nucleon-nucleon scattering requires in the ^5o channel that at leading order (LO) 
a constant contact interaction and one-pion exchange (OPE) be fully iterated. However, a large 
discrepancy exists between the resulting EFT LO \^, |6|, lZ| and partial-wave analysis (PWA) by, 
say, the Nijmegen group [0|, suggesting rather slow convergence of ChEFT expansion. One pos- 
sibility is that there may exist an unexpected infrared mass scale due to fine-tuning of quantum 
chromodynamcis (QCD) such that momentum (or energy) dependence must be somehow incor- 
porated into the LO short-range interaction, as opposed to WPC l|]. Reference [9| showed that 
such an infrared mass scale can be manifested by the inverse generalized effective range in the 
modified effective range expansion (ERE) for ^^o (to be defined model- independently in the pa- 
per), 2/r ~ 100 MeV, whereas its natural value would have been around the break-down scale 
of ChEFT — Mhi ~ rrifj ~ 600 MeV with ma being the mass of the a meson. This fine-tuning 
requires resummation of kr/2 to all orders, where k is the magnitude of the center-of-mass (CM) 
momentum. 

The large value of the ^Sq scattering length, defined by the regular ERE near threshold, is 
yet another, albeit much better-known, fine-tuning of QCD 10|]. The two fine-tunings do not 
seem to be correlated though, since r has more to do with the energy/momentum dependence of 
short-range forces while a is closely related to the constant part. Originally designed to deal with 
fine-tuning of a, the machinery of Ref. [l[, interestingly, facilitates resumming kr/2 to all orders: 
Introduce an auxiliary field, called dibaryon field and denoted by cf), that has the same set of 
quantum numbers as the ^5*0 partial wave — baryon number 2, parity even, spin 0, and isospin 1 — 
and the s-channel exchange of will bring the desired energy dependence to the LO potential. It 
is my goal to show in this paper how this machinery can be generalized to include systematically 
higher-order corrections, in particular, those of irreducible two-pion exchanges (TPEs). This is 



13[ to modify WPC and build a consistent and efficient power counting for 



part of our efforts [11 
chiral nuclear forces. 

This goes beyond those works that modify WPC in order to respect renormalization group 
(RG) invariance [gI. Il3l4l6l|. in which fine-tuning of momentum-dependent ^5o counterterms were 



not particularly considered. However, the findings of Refs. |l3l4l5(| make the fine-tuning of r appear 
less surprising than it would for WPC: In a natural chiral system the momentum dependence of 
^Sq counterterms would arise as 0((5/Mhi) correction to LO, as opposed to the underestimation 



of 0{Q'^/M^^ by WPC. 

Another line of investigation on applying the dibaryon fields to chiral nuclear forces can be 
found in Refs. 



17 



18| . in which the dibaryon fields are used for both ^^o and '^5i while only the 



spin-0 dibaryon is employed in my approach. In addition, the central Yukawa part of OPE [se e 



m. 



Eq. @] is treated as perturbation in these works along the line of so-called KSW counting 
whereas OPE is considered nonperturbative in this work [see the discussion below Eq. @.] 

For the notation to be more c omp atible with the literature, I depart from the convention 
adopted in our previous papers 



12 



13( 1 and denote the order of EFT amplitudes by their absolute 



size rather than relative size compared to LO. Therefore, nonperturbative LO will be labeled 
in the paper as 0{Q~^)^ which is the scaling for any nonperturbative, nonrelativistic scattering 
amplitudes, 0{Q/M\i\) corrections to LO as 0{Q^), and so on. 

Note that the WPC LO of "^Sq has actually another issue, which is not directly related to the 
aforementioned slow convergence: WPC fails to prescribe a quark-mass dependent counterterm at 
LO which is, however, required by RG invariance 20|, |2l[. As a by-product of the technique dis- 



cussed in this paper to tackle the slow convergence, the quark-mass issue will be solved altogether. 
I review in Sec.|lT]the theory without fine-tuning and introduce the unnaturally large, generalized 
effective range. I then show in Sec. IIIII how this fine-tuning can be incorporated by utilizing the 
dibaryon field, and then demonstrate the corresponding power counting up to 0[Q^) where the 
leading TPE needs to be accounted for. Finally a summary is offered in Sec. IIV[ 

II. ISSUES AT LEADING ORDER 



To motivate the employment of the dibaryon field, I briefly review the original theory that does 
not include it. Since most of points to be shown here were already made in the literature, the main 
function of this section is to establish the notation. 

The leading Lagrangian terms concerning the ^ ^o channel of NN scattering are [3, Il9l . 

72 
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C-NN = \{d^.^? - \ml-K'' + iV^ A^o + ^-^ N - I±NWaa,id,7ra)N 
2 2 V 2mNj 2/^ 



do{N^PaN)^N^PaN - Cr-^ 



1 - ^V^P. 



m^ 



{N^PaN)^N^PaN 



(1) 



+^ 



{N'^ PaN)'< N'^ Pa{^ - ^fN + 



h.c. 



where ttItt = 138 MeV, gA = 1.26, /tt = 92.4 MeV, niN = 939 MeV, and Pa is the spin-isospin 



projector for the ^5*0 channel: 

Pa = —r=T2Ta(y2 ■ (2) 

The quark-mass term proportional to Cq"^ is written in the so-called stereographic coordinates 



for TT [23|. Not only does it bring m^ dependence to the contact interactions but produces a 
non-derivative inrNNNN vertex. 

The leading ^Sq amplitude by WPC is resummation of OPE and a constant ^5*0 counterterm 
to all orders: 

V(~^)(q) = VY{q) + Co, (3) 

with the Yukawa potential and Co defined as 

W(.)^-i^^, Co.Co + ^^. (4) 

Here q = p' — p, with p' (p) being the outgoing (incoming) momentum in the CM frame, and 
a~^ = IQtt f^ / g\mj\f ~ 290MeV. If a~^ is chosen to be an ultraviolet (UV) mass scale, the 
resulting power counting is KSW. This is especially plausible if the degrees of freedom of the 
delta-isobar are integrated out, and hence the delta-nucleon mass splitting ~ 300 MeV becomes 
the break-down scale. However, I choose to work in a more general scenario that keeps open the 
possibility of incorporating the delta-isobar, that is, a~^ is considered an infrared mass scale in 
the paper: a~^ ~ Mi^. In addition, this allows for exploring situations where m^ becomes so large 
that mT^aTj ~ 1. 

The LO ^^0 amplitude can be written as 

T^ >{p ,p;k) =Ty{p ,p;k)+ — , (5) 

-L/L-o - Jfe 

where Ty is the fully nonperturbative iteration of Vy and 

X{P,k}-l+ (2vr)3 £; _ JL + ie ' ^^^ 



{2n)^E-^ + 



rriM 
U2 



le 



with E = k /rriN being the CM energy. Shown diagrammatically in Fig. [H the expansions of 
Xk = x{k; k) and I^ in powers of Vy suggest that Xk is finite while I^ has divergences oc ^^[/JqA + 
(iia-j^m^ ln(A/;Li)], where A is the UV momentum cutoff, ^ is an infrared renormalization scale, and 
/3o 1 are coefficients depending on the form of the regulator. 








x(p; k) = 
Ik = 



FIG. 1: Diagrammatic representation of x(p; k) and Ik- Here the solid (dashed) hnes represents the nucleon 
(pion) propagator, and the crossed circles represent no interactions. 

The divergent m^ dependence of Ik immediately calls WPC int o q uestion: Chiral-invariant 



1/Co cannot absorb such a chiral-symmetry breaking divergence [201. l2l|. (This is not entirely an 
academic issue that concerns only extrapolating lattice QCD results to the physical point of 771,^; the 
TTi^ divergence is related through chiral symmetry to one of the non-derivative irnNNNN vertexes, 
which may have phenomenological impacts.) Besides the obvious workaround that is the KSW 

2l| that promoting to LO the operator proportional to d2m'^ — 



counting, it was proposed in Ref. 



replacing I/Cq with l/(Co + d2rn^) in Eq. ([5]) — will be able to renormalize the m^ divergence of 
Ik- This is quite a striking statement because it is not immediately clear how l/(Co + d2rn^), as 
a fractional function, can absorb a linear divergence oc m^. As we will see, the dibaryon field can 
resolve this issue in a more transparent way by allowing its mass to be renormalized by iterations 

of Fy 

Although important, the short-range structure related to quark masses does not change the 
form of the ^^o amplitude as a function of the CM momentum, and hence does not help resolve 
the other issue at LO which is more phenomenologically urgent, namely, the slow convergence of 
EFT expansion in comparison to PWA. Thus, for simplicity in the qualitative discussion of the 
convergence issue, I allow for the moment Cq to have nontrivial m^ dependence before we proceed 
to serious calculations. With this caveat, one can use I/Cq to cancel both divergences of /^ and 
writes on-shell T^^ as 

where I^ilA is the /x dependent, finite part of I^ uM- After fitting to the PWA, one immediately 
observes the slow convergence, through the rather large discrepancy between T~^ and the PWA 
at low energies (see, e.g.. Fig. 2 of Ref. [13|].) For instance, the LO EFT predicts ~ 65 degrees for 
k ~ r7T,7r while the Nijmegen PWA gives ~ 40 degrees. 

However, the next-to-leading-order (NLO) EFT curve fits the PWA well, which means that 
the NLO amplitude, T^^\ is unexpectedly enhanced. Before any fine-tuning is considered, the 



RG- invariant chiral theory for ^Sq requires that (renormahzed) contact couphng constants Cq and 
C2 scale as fohows 13h15|. 119| , 

Att 1 „f? 47r 1 






and hence the ©(Q/Mhi) correction to LO ^ is 



Ci 



ITIN M^M^ 



C^ 



k'xk' 



v2 • 



(9) 



(10) 



To see more easily how Cq and C2 are linked to the phase shifts, I recast T^^^' + T^^' in the form 
of the modified ERE, 

47r 



7^(-i) _^ 7^(0) ^ 21, 



Y 



xl 



"^^ -SR + 1^' + ^4^(^) 



1 + 



A: 



(11) 



where 



47r 1 



a(/i) mAT C^ 
f 47r C2^ 



Mio, 



1 



(12) 
(13) 



2 m^ (Co«)2 Mhi • 

Unlike a(^), which depends on the renormalization scale, r is well defined and its value can be 
extracted from the phase shifts [9|: 

1 



r 
2 



1.55 fm 



(14) 



127 MeV 

The fact that 2/r <^ Mhi — 600 MeV signals that the data do not faithfully support the proposed 
scaling for r/2 in Eq. Q — the only avenue through which Mhi could have suppressed T^^' . To 
accommodate the empirical fact r/2 ~ 1/Mio, we need to develop a new scheme in which the ratio 
kr/2 are resummed to all orders such that rk^ /2 occurs in the new scheme at LO rather than 
appears as a subleading correction. 

The regular ^Sq scattering length, defined as the zero-energy value of the amplitude, is also 
unnaturally large: 015^ ~ — 24fm, compared with its would-be natural value of 0(1) fm. It is 
tempting to consolidate the two fine-tunings and to argue that they come from the same source, 
but ChEFT will not be able to do this because at fc = where aig^ is defined, rk'^ /2 does not 
contribute. In other words, as far as ChEFT is concerned, dialing r does not seem to have any 



effects on ai 



So- 



Promoting nothing but Co to LO seems an obvious way to achieve the resummation of rk/2. 



This was proposed in Refs. 



20 



21 



24| and was numerically shown to work well for a range of 



cutoffs 21[. By arranging some sophisticated runnings for Cq and C2, Ref. 25|] claimed that one 



26[ argued that iterating both Cq 



can obtain analytically a renormalized LO amplitude. But Ref. 

and C2 with different regularization schemes will lead to different results and that with a cutoff 

regulator the effective range allowed by the theory cannot be freely chosen — the so-called Wigner's 

bound 



23, 



28l |. So it is still unclear whether one can promote C2 alone without sacrificing RG 



mvariance. 

Besides the aforementioned technical difficulty in renormalization, it will be quite surprising if 
fine-tuning of C2 can be isolated without contaminating operators with four more derivatives. This 
is because multiple insertions of lower-order counterterms will generally renormalize higher-order 
ones through loops. Before fine-tuning being considered, C2n — the coefficient of the four-nucleon 
operator with 2n derivatives — are loosely correlated by Mio and Mhi, through dimensionless coef- 
ficients, 62n, that are undetermined but are of 0(1) (l3[: 



%(/"+/") 



Air 



72n 



n 



ruN MI'+^Ml'. 



{p'^'+p 



/'in-. 



(15) 



lo 



^hi 



(Note that Ref. [IJ] proposed different scalings for C2n, hence a different type of correlation ensues.) 
62 being tuned towards larger values while ^0 remaining fixed effectively lowers Mhi, and such change 
of Mhi will propagate to C2n with n ^ 2. 

Now I recklessly assume that the above correlation is still pertinent even when C2n{n ^ 1) are 
tuned to be so large that Mhi is replaced by M\o in Eq. (jlSp . It then becomes apparent that all 
of C2n will be equally important. But an EFT with infinitely many unknown parameters at LO is 
not meaningful, unless we impose a stronger correlation among 92n such that C2n are determined 
by a finite number of LO parameters. Consider the following correlation of C2n by two parameters 
at LO, Oo/Mio and 6'2/Mio, 

47r Oq ^ Att f 00 Y^^ M2 ^" 



Co 



C2r, 



At tree level p' = p = k ^-^ Mi^, the sum of all C2n^^" is 
^' Att 1 



A/io 



5]C2„fc2« 



n=0 



THN Mio/^0 - 02kyMy, 



l + O 



k' 



MhiA/i 



lo 



The above summation resembles a tree-level s-channel exchange of in NN scattering: 

ay'^ Att 1 



with 



^ + A 



47rA 
mAry2 



rriN (jJ2L 



fc2 



Mlo: 



A _| 



Att 



2 9 



Af, 



lo 



(16) 



(17) 



(18) 



(19) 



Here a = ±1, A is the mass splitting between cf) and two free nucleons, and y is the (pNN coupUng. 
As first shown in Ref. pj and to be reiterated in the next section, the s-channel (f) exchange will 
bring about the desired LO amplitude in the form of Eq. (jlip. which is the ultimate justification 
for me to have chosen such a correlation as Eq. (|16|) . 



III. POWER COUNTING WITH THE ^Sq DIBARYON FIELD 



The first few Lagrangian terms involving cf) are 



m 



Cf = ac/.' ■ I iSlo + j^ + A I lA + y ( 4>iN' P„Af + h.c. 



+ * "' (rTl#) *'■'*+'"= "'- (iT^) (^i""^"" + '"■) ""' 



where ^^ is the covariant derivative for an isovector field: 

2 \-l 



= a,0-(i + ^) (^xl^)xc/,. (21) 



4/2y \u " 2A 

Here I have normalized cf) so that cr = ±1. The (i2 term are chiral-symmetry breaking and, as will 
be seen, it needs to be at LO because of the fully iterated Yukawa potential. Terms that do not 
explicitly involve cf) but will be needed in the paper are already shown in Eq. ([T|). 

A. O(Q-i) 

Following the argument that leads to Eqs. (J17p and (jlSp . I revise WPC such that the LO 
"short-range" potential is represented by an s-channel exchange of cf): 

y(-i) = V^{E) + Vy , (22) 

where 

VJE) = ^ ^ . (23) 

Having both mass and kinetic terms of cf) at LO means that I have chosen V(jy{E) ~ Vy ~ 
-^iii—r-7-TT77\, which means that for renormalized A, do, and v: 

A«~dfm2~^ and [y^f ^ ^^ . (24) 

niN iTij\f rriN 



The new LO potential is computationally equivalent to making Cq of the dibaryon-less theory 
energy-dependent. The new LO ^Sq amplitude is then found to be 

T^^^) =Ty+ ^^ 



a^±^ + a^^' 



Ty ^^ ^' 



y^m^ (25) 



with the generalized scattering length and generalized effective range defined for A — )• oo as 

1 47r f CT A-^ 1 cfd^ \ 

+ ml-^ ^ ~ Mio , (26) 



a(^) TUN \{y^Y ''(y^^^ 
r Att a 1 



where 



2 TTlTV fT^NV"^ -^^lo 



A'^ A .^ . rf^ -J2 «„, , /A\ ..,., .« 



(27) 



^ (P?j2 - ^^ " ^0^ ' "" (7?)2 = ""^ " ^1"- ^'' v;^J ' ""^"^ y =y- (28) 

Thus, we arrive at the desired scalings for a and r. For a finite value of A, r has residual A 
dependence that vanishes at the rate of 1/A. The details of the numerical calculations will be set 
up later, I only remark that for A = 800 MeV, I find r/2 = 1/115 MeV, consistent with its value 
stated in Ref. [9|. 

We also see the necessity of having 6,2^1^ at LO; its assignment is to subtract the jtt,^ divergence 
of Ik- As a consequence of defying naive dimensional analysis, ^2"*-^ gives rise to unsuppressed, 
non-derivative Tnr<j)(j) coupling: 

71-2 

^n-Kip^, = -d2ml—^(t)^ • (f) . (29) 

The form of the LO ^Sq wave function in coordinate space will be needed for later use. Since they 
are somewhat out of the main line of the physics, I relegate the technical details of its construction 
to Appendix 1X1 

B. O(Q0) 

The most general dibaryon Lagrangian is bound to have many redundant terms because there 
will not be enough observables to pin them down, due to the fact that does not correspond to any 
particle appearing in asymptotic states. I choose to minimize the number of 0-related operators 
and to have four-nucleon contact operators be responsible for improving short-range interactions 



at subleading orders. This choice means that after the exchange taking away the dominant part 
of short-range interactions, Cq of the dibaryon Lagrangian represents higher-order effects: 

Co-—^- (30) 

Note that I have shghtly modified the scheme of Ref. []J in which Co was put on equal footing with 
the cj) exchange. 

It is convenient to expand formally bare low-energy constants (LECs) to reflect the fact that 
even though the number of physical inputs must stay the same, their RG running may change at 
each order, 

A^ = A(-i) + AW + A(i) + • • • , (31) 

df = 4"^) + df + d« + • • • , (32) 

yB^yM)+y(o)+y(i) + ...^ (33) 

C^ = cS'^ + ci,'^ + ---, (34) 



where the expansions are in powers of 1/Afhi. For each parameter, the superscript of the leading 
term in its expansion marks the order it starts to contribute. For instance, since Cq is the first 
term in Eq. p4p . 0{Q^) will be the order Co occurs for the first time. However, in order to improve 
the readability of the manuscript, I make a few exceptions and will drop the superscript ^^^' for 
the first term of A, d2, and y. 

NLO potential V^^' consists of only contact interactions and corrections to A, d2, and y: 

y(0) = C7(°) + 2 (y (0) + 4°)m2 ) Y± - a (a(o) + df^ml + df m^) [^ V . (35) 

J \ y / 

Despite the energy dependence of V^, the technique shown in Appendix B of Ref. [13(] is still useful 

in evaluating insertions of V^^' . One can find a single insertion of V^^' to give rise to the generalized 

shape parameter, in addition to m^ and m^ corrections to 1/a and r/2, respectively: 



r(o) = ^ 



;i)(o) + M(°^2+~^^4 



xl 



rriN ( 1 I 7,2 47r jR^"^ 



(36) 



10 



where 

a) rriNy^ y y-^ V J 

^2 = ^Q^-. (39) 

Here I have chosen A'"-*, dg , and y^^' to be such that the chiral invariant parts of 1/a and r/2 
and m^ part of 1/a remain their LO values: 



y2 I y2 ^ ^ 



2<tA 

y2 



^^^^ + 2— -^1=0, (40) 

0, (41) 



y2 I y y Ai A 



,^f.^,<^U„. (42) 



r*- ^^ + T^^' can be rewritten in the form of the modified ERE: 



T(-i) + tW = Ty - i^ , ,,. ^j,. ,_. + O f #T("i)') . (43) 



m^_i + |fc2+^^fc4 + J^ji? VA4i 



Power counting (j30p is then equivalent to estimating V2 as 



Mhi 4 M2 Mhi ' 
which is compatible with the value extracted in Ref. 



^ 1 f^ 1 

^2 - T7-ir - .^2 ^, ' (44) 



^^=5ifeV' "'*^ i=127MiV' (^^) 



and the value by this work for A = 800 MeV, 

r^/4 r 1 

^2 = , with - = . (46) 

693MeV ' 2 llSMeV ^ ' 

c. o(gi) 

1. Residual counterterms 

It is, if only academically, interesting to ask how counterterms will scale if TPEs and higher-order 
multiple-pion exchanges are completely turned off while the strength of OPE remains unchanged, 
which can be achieved by taking 1//^ — )■ but keeping niN/f^ fixed. Higher-order counterterms 

11 



in this scenario, referred to as "residual counterterms" in Ref. [Ij], will be responsible for all the 
subleading corrections, and the modified ERE is expected to be valid to all orders and to acquire 
j^2n ^gj~j^g beyond V2k'^, 

rp _ ^"^ Xk f'471 

n=2 

and I wish to find out how Vn (for n ^ 3) will scale. 

Next-to-next-to-leading-order (NNLO) amplitude T'^' includes two insertions of Co, which can 
be compared with a single insertion of C2(p^ + p' )/2: 

^.(1) _i^ X| (~j^6 V2k^ \ 

y a 2 misi k J 

with 

and • • • referring to m^ corrections to l/o, r/2, and V2, which are not pertinent to the present 
discussion. Cg is not running with A, but it is nonetheless renormalized by a term quadratic 
in Cg [Cq is related to V2 through Eq. (|39p .] With fine-tuning having been accounted for by 
the resummation at LO, it is reasonable to expect naturalness to retain its power in counting. 
Therefore, the two terms contributing to ^3 must have similar sizes, resulting in 

_ 1 47r 1 

vh ~ — 5 — ^ and Co ~ tt ■ (50) 

More generally, one can show that Vn+i will have contributions, among others, from n insertions 
of Co and one insertion of C2n-2{p^" +p' ")/2: 

..+i = ^ + ^C,„_,- + ..., (51) 



which leads to 



1 „ 47r 



^"+1 ~ .fU+U.n ' ^2n-2 ~ ,.„-!,.„, " (^2) 



Ko ^U "^^ Mo" ^hi 



^. Two-pion exchange 



Let us turn to the leading TPE, V27r. Throughout our efforts 



11 



13( 1 to modify WPC, we have 



taken the position that the standard chiral counting does not need to change for pion-exchange 



12 



diagrams [4], which essentiahy describe long-range physics through non-analytic functions of mo- 
menta. This means that since V2tt is suppressed by C'[Q^/(47r/^)^] relative to OPE, the single 
insertion of V2tt into the resummed LO amplitude remains an C'[Q^/(47r/7r)^] suppression.^ 

But one still needs to to determine what counterterms are needed to renormalize the UV part 
of the single insertion of V2tt- They are called "distorted- wave counterterms" in Ref. [13|] because a 
single insertion of V27r is equivalent to the matrix element of V27r between the LO wave functions — 
the distorted wave for the LO potential [ij, [iGl] — 

(V'fc|^2.|V'fc>=47r / drr^i;l{r)VJ^\r). (53) 

Near the origin, ipkii") can be expanded in powers of (fcr)^ [see Eq. (|A13p ] and V2tt oc 1/r^; therefore, 
one can find the superficial divergence of {ipk\y2n\ipk) to be 

An [ drr''^l{r)vi^{r) 

2 (54) 



oc 



^^ * {pqM^ + pi7k^K-^ + p2r'^k^A^ + P3?'^k^ In A) + F.T. 



V7^ - I, 



where "F.T." refers to finite terms and pn have at most logarithmic dependence on A. For each 
/c^" only the highest power of A is shown. 

To identify the needed counterterms, we first notice that C2(p^ +p' )/2 produces k^xt/i^^ " 
Ik)^, as suggested by Eqs. ([18|) and (jl9]) . With a little algebra, one can show that other divergences 
with lower powers of k'^ can be subtracted by corrections to A, y, and Cq. With every piece put 
together, V^^' has the form: 



V^'^ = V2n + ^iP^+ P') + Ao + XiV^iE) + X2V^{E) + X3V^{E) . (55) 



It would be such an unnecessary bore to write down the expressions of Aj in terms of A'-^^ y , 
etc., because in practice all we need to know is that Aj are independent of energy or momenta. 
I have also stopped pursuing the complete understanding of m^ dependence, which will be dealt 
with in future publications. 

It simplifies tremendously the work of establishing power counting that V27r does not demand 
for renormalization purpose more counterterms than the residual counting ()52p provides. I ex- 
trapolate this observation to any irreducible multiple pion exchanges and conclude that all of ^Sq 
counterterms are prescribed by power counting (I52p . 



However, I would like to remind the reader that a different point of view towards numerical factors of tt in chiral 
counting can be found in Ref. [30|. 

13 



D. Results 

Although the formal expressions shown in previous subsections reveal the renormalization and 
the analytic structure of the amplitude at each order, a complete analytical calculation is still un- 
likely for quantities like x(pj k) and Ik cannot be computed analytically. So numerical calculations 
are necessary. 

The potentials are regularized with a separable momentum-space regulator: 

Va{p',p; E) = exp L^l-±^ j V{p',p; E) . (56) 

While the LO amplitude is generated nonperturbatively by solving the Lippmann-Schwinger equa- 
tion for the LO potential ()35p . 

2-(-l) ^ y(-l) + y(-l) Q rp{-l) ^ (57) 

where G is the Schrodinger propagator, the subleading corrections are calculated through pertur- 
bative insertions of y(°) ([35]) and F^^^ ([55]), 



rw = (i + r(-i) G) y(i) (g t^-^) + 1 j , (58) 

^(1) = (i + r(-^) g) [y(2) + yd) (g + gt^-'^ g) v^'A (gt^-'^ + 1) , (59) 



in which I will adopt from Ref. [J] the delta- less version for V2n- 

Figure [2] shows the EFT results for ^5*0 phase shifts up to and including 0{Q^). In the left 
panel, the LO is obtained by fitting to the PWA points at Tiab = 5 and 15 MeV. At 0{Q^) and 
0{Q^), 25 and 50 MeV are added, respectively. The bands are generated by A = 0.6 — 2 GeV. 

. er_than one tenth degree at, e.g., Tiab = 130 MeV. 
, the new formulation fits much better to the 
PWA and converges rapidly, at the cost of one more short-range parameter at each order. The 
breakdown of convergence around Tiab — 200 MeV is not surprising, for the delta-isobar is not 
explicitly considered. 

The small correction provided by the leading TPE reassures its perturbative nature. This is in 



Above A = 2 GeV, the cutoff variation is smal 
Compared with the dibaryon-less theory |13l . 



contrast to the WPC-based study of Ref. 



3l|, which suggested that in order to have perturbative 



32| 



multiple-pion exchanges the cutoff needs to take smaller values, and the speculation of Ref. 
that the delta-isobar may be the source for the slow convergence. 

To get an idea of how much the fit can be further improved, I fit to PWA points at higher 
energies, and the results are shown in the right panel of Fig. [2j There the fitted PWA points are 
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50 100 150 

Tiab(MeV) 



50 100 150 
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200 



FIG. 2: (Color online) ^5*0 phase shifts as a function of laboratory energy. The black dots are from the 
Nijmegcn PWA 8|. In the left panel, the light blue, green, and dark orange bands are 0{Q~^), 0{Q^), and 
0(Qi) calculated with A = 0.6-2 GeV. In the right panel, ©(Q") (green dot-dashed hue) and 0{Q'^) (dark 
orange solid line) are plotted with A = 1 GeV. 

at Tiab = 50 MeV for 0(^0) and Ti^b = 150 MeV for 0{Q^). Since the cutoff dependence is no 
longer a concern, only A = 1 GeV is used for 0{Q^) and 0{Q^). 



IV. SUMMARY 

I have considered the EFT expansion for ^ Sq of A^A^ scattering in which the generalized effective 
range r/2 is counted as an independent, infrared length scale due to fine-tuning of the underlying 
theory. The new expansion was made possible by an auxiliary, dibaryon field that has the same 
quantum numbers as the ^5*0 partial wave j20( |. At LO, the fine-tuning in question is incorporated 
by iterating the s-channel exchange of cf) to all orders. A by-product of employing the dibaryon 
field is the capability to easily cancel the m^ divergence that stems from the fully iterated OPE. 

I chose to minimize the number of 0-dependent operators and to use four-nucleon counterterms 
to account for subleading short-range forces. If irreducible multiple-pion exchanges were hypothet- 
ically turned off, power counting of the so-called residual counterterms could be considered. C2n 
would appear in the 0{Q"'^^/M^^ ) corrections to LO and would scale as 

An 1 



O 



2n 



mj, M^M'^, 



n+l ■ 



(60) 



This counting was actually found to provide enough counterterms to absorb the divergences of 
TPEs, when they are turned back on and are inserted to the LO amplitude. Therefore, the 
above power counting is the final answer we were looking for. The numerical results showed much 
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improved convergence of the EFT ^Sq phase shifts. 
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Appendix A: LO wave function 

I follow the technique developed in Ref. [201 ] to obtain the short-distance behavior of the LO ^5o 
wave function. In the limit A — )• oo, the LO potential has the following formal coordinate-space 
form: 



y(0)(f) = y4^)5(3)(f) + FY(r), 

where Vyi'r) is the Fourier transform of the Yukawa potential and V^ is defined in Eq. 
in-state, 5-wave wave function formally satisfies 
1/^2 2d 



+ Vy-E 



Mr) = -V4E)iJk{0)S^^Hn , 



(Al) 
The 

(A2) 



mj\f \ dr"^ r dr , 

and can be written as a linear combination of the regular and irregular solutions to the Schrodinger 
equation for Vy, 



Mr) = aik)Mr) + b{k)nk{r) , 



where Jk{r) and T-Lk{r) are normalized so that they satisfy 

1 



d'^ 2d,,, 
rriN \ dr'^ r dr 



1 



d^ ^ d , ^^ 
ruN V dr^ r dr 



Jk{r) = 0, 



(A3) 

(A4) 
(A5) 



While Jk{r) — ^ joi^r) for r — t- 0, where Jq{x) is the zeroth spherical Bessel function, "H^ has the 
following form for r — >■ 0: 



n,(r) 



A-K 



B{kr, K-j^r) — 2KTrA{kr, n^^r) In (//r) 



(A6) 
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where Kj^ = w,^a7r and A{x, y) and B{x, y) are dimensionless functions that are analytic at x, y = 0. 
Using the above expression one can obtain the expansion oiTikir) in powers of (kr)'^ and/or (K.^r)'^ 
near the origin. 

For any cutoff regulator, the delta potential gets smeared away from the origin, up to a dis- 
tance characterized by 7^ = A~^. Solution ()A3p in fact governs the "outside region", r > IZ. 
Reference |20( | showed that the singularity of T-LkCR-) for 7^ — t- can be related to the diver- 
gences of /fc. To see this, notice that Eqs. ()A4p and (jASP indicate that certain linear combination 
of Tikir) and Jk{i~) makes up the S-wave interacting Green function for the Yukawa potential: 
Gyir; E) = (r, yQ{0, (^)| {E — Hq — Vy + ie)~ \x = 0) with Hq being the free two-nucleon Hamil- 
tonian and Yq{9,(J)) the S-wave spherical harmonic. The divergences of Gy(0;£') is completely 
described by T-LkiJ^ — ^ 0) and does not depend on the Jk part, since Jk{r) behaves well near r = 0. 
On the other hand, Ik is precisely Gy(0; E), most easily seen from its diagrammatic representation 
in Fig. [TJ Now we can identify the divergences of Ik, — /3oA — /3iK7r ln(A//_f), with TikiTZ — ?• 0). 
This means that following subtractions in Eqs. (j26p and (j27p . one can also use Vr to subtract the 
singularity of 7^fc(7^ — ;> 0): 



Substituting Eq. (|A3P in Eq. (|A2p and applying Eqs. (|A4p and (jASP on the left-hand side of 
Eq. (|A2l), we find 

b = -V^{E)i;km =b[l + nk{0)] . (A8) 

On the other hand, V'fc(O) is generally related to the LO off-shell T-matrix by 

Using Eq. ([5]), with Cq being replaced with V0(E'), in the above equation, one finds 

^.(0) = ^^ = — ^^^^^V^ . (AlO) 



(All) 



Using the above expression and Eq. (IA7p in Eq. (|A8p gives 

, Xk 



vr' - h 



k 



a 
5 



V+^^(0)]=-^(4 + l^')- (^12) 
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It is worth noting the difference between this LO wave function and that of the dibaryon-less 



theory: while b stays t 
dibaryon-less theory [l 



le same, a/b now has k dependence, in contrast to being a constant in the 
■ Il3l|. a consequence stenmiing from the fact that the LO contact interaction 



now provides two inputs for the outside wave function. 

We can now put these back into Eq. (|A3p to have a more clear expression of the short-distance 
behavior of ipk (^) ■ 



Mr) = --^^^[lA{kr,.^r) 



2K^ln(/ir) - - -I- -/c^ 
a 2 



B{kr, K^r) } . (A13) 



This expression is in principle accurate only for TZ <^ r <^ k^^, and the specification of the 
regularization scheme is expected to change the details of the wave function near r ^ IZ. But 
this does not invalidate the qualitative statement I made in Eq. (I54p regarding the divergences of 

('0fc|^27r|'0fc>- 
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